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The unsurpassed sensitivity and resolution of the Square Kilometer Array (SKA) will make it possible for the 
first time to probe the continuum emission of normal star forming galaxies out to the edges of the universe. This 
opens the possibility for routinely using the radio continuum emission from galaxies for cosmological research as 
it offers an independent probe of the evolution of the star formation density in the universe. In addition it offers 
the possibility to detect the first star forming objects and massive black holes. 

In deep surveys SKA will be able to detect Hi in emission out to redshifts of z ~ 2.5 and hence be able to 
trace the conversion of gas into stars over an era where considerable evolution is taking place. Such surveys will 
be able to uniquely determine the respective importance of merging and accreting gas flows for galaxy formation 
over this redshift range (i.e. out to when the universe was only one third its present age). It is obvious that only 
SKA will able to see literally where and how gas is turned into stars. 

These and other aspects of SKA imaging of galaxies will be discussed. 



1. Introduction 

More and more evidence supports the paradigm 
that galaxies form hierarchicallly, i.e. that 
present day galaxies have been built up over cos- 
mological time scales by the coalescence of many 
smaller, less massive galactic systems {415716518^1 
11911) . For a particular galaxy this process is de- 
scribed by the so-called merger tree: the accumu- 
lation of dark matter mass over time into a final 
halo. These merger trees are extracted from nu- 
merical comological simulations and only describe 
the merging of dark matter haloes in a robust 
manner. 

Many properties of galaxy populations includ- 
ing current star-formation rates, luminosity func- 
tions, the Tully-Fisher relation, morphology seg- 
regation, and the variation of star-formation his- 
tories with environment are explained by semi- 
analytical, hierarchical models. However, the role 
and fate of the gas and stars in this merging pro- 



cess is not well understood due to an incomplete 
understanding of the relevant physics, the sim- 
plicity of the hydrodynamical codes and the lack 
of numerical resolution. The behaviour of the 
baryons in these merging dark matter haloes is 
often described with semi-analytical techniques 
|5:il7ll9l2T)j (see also Baugh et al. (this volume)). 
In this scheme, the physics of galaxy formation 
is encoded in a set of rules that are either phys- 
ically motivated or have been designed to repro- 
duce the results of numerical simulations e.g. the 
consequences of galaxy mergers. A merger tree 
(the "family history" ) of a dark matter halo is ei- 
ther generated using Monte Carlo techniques or 
is extracted from an N-body simulation. The up- 
permost branches of the tree are filled with gas 
that is shock heated to the virial temperature 
of these progenitor halos, and a set of rules is 
followed to describe the subsequent evolution of 
this gas. Apart from the issue of properly mod- 
elling the fate of the baryons, it is also true that 
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the typical shape of a galaxy merger tree has not 
been confirmed robustly from observations. Ob- 
servationally, the galaxy mass functions and the 
merger rates as a function of redshift are largely 
unknown. 

In this picture there still are several key ques- 
tions that need to be addressed. These are: 

- When and how quickly did the gas settle 
into the dark matter potential wells? 

- Where, when and how did the gas get con- 
verted into stars? 

- When and how did the merging process of 
small galaxies develop? 

- What is the relative importance of gas ac- 
cretion versus galaxy merging? 

- Are galaxies still (trans)forming today? 

In order to address these questions it is necessary 
to have a census of galaxies and their properties 
over a large range of environments and as large a 
possible range of cosmic time or redshift. While 
the stellar content of galaxies can be probed using 
optical and NIR imaging and spectroscopy using 
the best facilities on the ground and in space, the 
only way to probe the cold gas content of galax- 
ies is to use a radio telescope, either at very short 
wavelengths to measure molecular line emission 
(see also Blain et al. this volume on molecu- 
lar lines), or at long wavelengths to measure the 
Hi emission using the 21-cm line. To do this out 
to redshifts of z = 3, since when considerable evo- 
lution may have taken place (83137) . requires a 
radio telescope tenfold larger (in collecting area) 
than any of the existing facilities. The SKA is 
such a telescope and the purpose of this chapter 
is to review the current state of affairs, from both 
the observational and theoretical points of view, 
and to outline how the SKA will be instrumental 
in answering the above questions to which we do 
not yet have firm answers. 

Equally as important, SKA can probe the star 
formation rates in galaxies using the radio con- 
tinuum emission. The radio continuum emission 
provides an estimate of the star formation density 
with redshift independent of other tracers such as 



optical or ultraviolet light which is affected by ex- 
tinction. The caveat is that the relation between 
the non-thermal radio emission and the star for- 
mation rate is assumed not to vary significantly 
with cosmic time. 

The framework for the discussion in this chap- 
ter will be the hierarchical structure formation 
models of Baugh et al. (this volume), who use 
the technique of semi-analytic modelling to pre- 
dict the evolution of the HI mass function and 
the number of star forming radio galaxies, both 
signposts of evolution to be probed with SKA. 
This follows a widely accepted idea of how galax- 
ies form, based on the success of the Cold Dark 
Matter (ACDM) simulations. 

Although ACDM has proven to be very suc- 
cesful in reproducing the WMAP results and re- 
sulting large scale structure, SKA will observe 
the baryonic component of the universe directly. 
Predictions of what we should observe are still 
wide open, as recent hydrodynamical simulations 
indicate that some of the basic assumptions of 
the semi-analytic models, i.e. that the gas gets 
heated quickly to the virial temperature, may not 
be correct. In the semi-analytic modelling the as- 
sumption is that the gas gets immediately shock 
heated to the virial temperature as it flows into 
the dark matter potential wells. Recently Bin- 
ney ^01, Katz et al. [2] and Murali at al. [HE] 
have pointed out that hydrodynamical simula- 
tions suggest that most of the gas is not heated 
to the virial temperature but accretes slowly in a 
cool phase from the surrounding large scale struc- 
ture filaments. The consequence is that galaxies 
grow through accretion of gas rather than merg- 
ing of dark matter halos. At redshifts beyond 
z ~ 2 this cool mode of accretion appears to dom- 
inate over merging as the main process for galaxy 
formation. Because these two scenarios differ in 
the evolutionary scenario of the gas, SKA is the 
ultimate telescope for distinguishing these from 
one another. 

We will first discuss the properties of radio con- 
tinuum emission from galaxies and in this context 
the first luminous objects in the universe. Then 
we will discuss the subsequent evolution of galax- 
ies following the theoretical scenarios, and finally 
connect this with the evolution of galaxies which 
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is still ongoing and observable at the present time 
in the local universe. 



2. The radio continuum emission from 
star forming galaxies 

When and how the first episodes of star forma- 
tion took place remains a key question of mod- 
ern cosmology. Much effort has been devoted 
to measuring the integrated cosmic star forma- 
tion rate (SFR) of the universe as a function of 
redshift or look-back time, initially at ultraviolet 
wavelengths and more recently in the mil- 
limetre and radio regimes |12I38| . Unlike ultra- 
violet or optical light, radio continuum emission 
is unaffected by dust, which eliminates the need 
for uncertain extinction corrections. SKA will 
detect the radio continuum emission from very 
large numbers of AGN and star-forming galax- 
ies at redshifts z > 3. However, we caution that 
the conversion between a measured radio contin- 
uum flux and a derived a global SFR requires the 
adoption of canonical scaling relations [22] which 
remain valid out to at least z ~1 |36I3) but, as we 
discuss below, may not be directly applicable to 
primordial star- forming systems at z > 3. 

In the models considered here, the radio contin- 
uum emission from star-forming galaxies is com- 
puted using the spectro-photometric code GRASIL 
written by Silva et al.[J§|- This code takes star- 
formation histories generated by GALFORM, includ- 
ing bursts of star formation, and produces a full 
spectral energy distribution, ranging from the far- 
UV to the radio. The radio emission includes 
both synchroton radiation by electrons acceler- 
ated in supernova remnants and free-free emis- 
sion from HII regions. The plots presented in 
this section are for a model very similar to that 
of Baugh et al. (this volume). Note that be- 
yond redshifts of z > 6 inverse compton losses 
against the cosmic microwave background will be- 
gin to affect the synchrotron emission. The non- 
relativistic, thermal emission is not affected [IT)] 
so that radio continuum emission still is a very 
good tracer of star formation at high redshifts, 
albeit in another regime. 
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Figure 1. Predicted 1.4 GHz flux densities for 
star-forming galaxies (with SFRs of 1, 10 and 
100M Q /yr) at different redshifts. The horizon- 
tal line at 100 nJy represents the nominal SKA 
detection limit for an 8 hr integration. 



2.1. The first luminous objects (z>5) 

Little is currently known about the properties 
of star-forming objects at redshifts above z^ 5 
but the next few years are likely to see rapid 
progress in this area. 

Several recent results suggest substantial star 
formation at high redshift |5] . Studies of the stel- 
lar 'fossil record' of nearby galaxies for example 
(e.g. 00]) imply that the star-formation rate in 
the most massive galaxies increases with redshift 
and peaks at z=3 or higher. The recent detection 
of QSOs at redshifts with z>6 |2H]j in which the 
gas appears to have solar or above-solar metal- 
licity |H5] suggests that the massive first stars 
in the Universe had already formed by z=8 or 
even earlier. Observations with the HST and 
ground-based optical/IR telescopes have begun 
to provide direct measurements of the number of 
star-forming galaxies at z ~ 5-6 [14198147181) . 
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These studies find typical surface densities of 
~ 100 deg~ 2 for star-forming galaxies at z ~ 6. 
So it appears that the global star-formation rate 
at z ~ 6 is roughly similar to that at z ~ 3- 
4, though with substantial uncertainities because 
the distribution of high-redshift galaxies appears 
to be clustered on scales comparable to the di- 
ameters of the individual fields studied to date 

M- 

At even higher redshift, the detections of indi- 
vidual (gravitationally-lensed) galaxies at z ~ 7 
and z = 10 [HE] have recently been reported, 
but the density of star-forming galaxies at z > 7 
is essentially unknown at present. 

SKA will be a sensitive probe of high-redshift 
star formation through the detection of radio 
continuum emission from star-forming regions 
and the radio afterglows of gamma-ray bursts 
(GRBs), many of which are now known to arise 
from the supernova explosions of massive stars. 

Models for dwarf galaxies forming in the vicin- 
ity of high-redshift QSOs imply likely star- 
formation rates of 30-60 M©/yr for these ob- 
jects |57) . Deep continuum surveys at 1.4 GHz 
with SKA will reach detection limits of 100 nJy 
in 8hr, and as low as lOnJy in longer integra- 
tions of — lOOOhr This should allow objects 
with star-formation rates as low as 10 M©/yr 
to be detected at redshifts up to z— 4 (8hr) and 
z— 11 (in 1000 hr) respectively (see Fig.^l. Indi- 
vidual globular clusters forming at high redshift 
may also be detectable by SKA, though there are 
currently few available predictions for the star- 
formation rates in such objects. 

2.2. What will SKA see? 

SKA will be able to detect very faint (10- 
100 nJy) radio continuum sources, i.e. 2-3 orders 
of magnitude fainter than those seen by present- 
day telescopes. It is therefore important to keep 
in mind that simulations of the likely surface den- 
sity, angular size and redshift distribution of these 
faint sources require extrapolations from our cur- 
rent knowledge (as well as assumptions about the 
cosmic evolution of each population). Neverthe- 
less, the work done so far (e.g. |4fil3f>l94| ) gives 
us some valuable ideas about what to expect. At 
100 nJy, the surface density of continuum sources 



will be very high (-200,000 deg" 2 ), with "nor- 
mal" , starburst and active galaxies all contribut- 
ing over a wide range in redshift. Simulations by 
Hopkins et al. imply that very high-redshift 
(z > 5) AGN and starburst galaxies will have sur- 
face densities of ~l,000deg -2 and — 7,000deg -2 
respectively in surveys with a detection limit of 
100 nJy at 1.4 GHz. 

2.3. Distinguishing starbursts from AGN 

Although star-forming galaxies dominate the 
radio source counts at flux densities below about 
1 mjy, active galaxies still account for a signif- 
icant fraction of the population at levels of a 
few n Jy |39| . The problem of distinguishing 
star-forming galaxies from AGN is discussed by 
Garrett et al. |35| . Several parameters mea- 
surable by the SKA itself (including brightness 
temperature, morphology, time variability and 
possibly polarization) can be used, and Garrett 
et al. I35| show that radio continuum observations 
with high spatial resolution (~ 25 mas) can reli- 
ably distinguish AGN from star forming galaxies 
at faint flux levels. This would require the SKA 
to have good sensitivity on baselines of at least 
3000 km. High spatial resolution is also required 
because of the likely source density at flux levels 
of — lOnJy, which may be as high as 3arcsec~ 2 

2.4. Probing high redshift star formation 
through GRBs 

Because of their transient nature, which makes 
them easy to recognize, and because they appear 
to be associated with the supernova explosions 
of young, very massive stars, GRBs may be one 
of the most direct tracers of star formation at 
very high redshifts. SKA, with its wide field of 
view, should be able to detect the radio after- 
glows of many 'orphan' GRB events in which the 
jet is misaligned with the observer's line of sight, 
and thus remains undetected in gamma-rays (e.g. 
Frail et al 29 ) . Bromm & Loeb \H>\ predict that 
at least 50% of GRBs originate at redshift z > 5, 
and Berger et al.[Hj find that GRB host galax- 
ies have typical star-formation rates of at least 
100M Q /yr. In the next decade, therefore, GRBs 
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are likely to become increasingly common probes 
of high rcdshift star formation up to (and per- 
haps even within) the epoch of reionization. 

2.5. Redshift measurements and estimates 

Although SKA will be able to detect many 
millions of star-forming galaxies over a very wide 
range in redshift, recognizing the highest-redshift 
objects will not be straightforward. For most 
star-forming galaxies with z<3, redshifts should 
be available from the large-area HI surveys which 
SKA itself will carry out. Alternatives at higher 
redshift include: 

(i) HI absorption lines see also Kanekar and 
Briggs, this volume). 

(ii) CO emission-line redshifts, using either the 
low order transitions observed with SKA, or 
higher-order transitions observed with ALMA, 
(see Blain et al. this volume) 

(iii) Radio 'photometric redshifts' based on the 
spectral slope between 1.4 GHz and 350 GHz |95] . 
This would require matching high-frequency ob- 
servations from ALMA. 

(iv) Photometric redshift estimates from wide- 
field multicolour optical/IR images, or spectro- 
scopic measurements in the optical/IR for small 
samples of particular interest. 



3. How different is high redshift star for- 
mation? 

In this section, we examine the assumptions 
underlying the use of the radio continuum to 
measure SFR, and assess their validity at the 
high redshifts where SKA will be able to make 
a unique contribution. We also outline new sim- 
ulations which will enable us to better interpret 
the planned deep surveys. 

Because of the far-infrared/radio (FIR-R) cor- 
relation (e.g. 02]), it is straightforward to show 
that the radio luminosity in evolved luminous 
star bursts varies only with SFR |22]. The radio 
continuum emission below about 10 GHz is domi- 
nated by non-thermal emission, which is thought 
to be proportional to the supernova rate i^snr- 
Although the mechanism linking ^snr to the non- 
thermal radio luminosity is complex, depending 
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Figure 2. 1.4 GHz radio luminosity (log) plot- 
ted against total far-infrared luminosity (taken 
from Bell 2003 [S]), showing the curvature at low 
luminosities. Normal star-forming galaxies are 
plotted as circles, intensely star-forming galaxies 
as stars (filled stars are ULIRGS, open stars are 
starbursts). Triangles are BCDs from a compari- 
son sample. The dotted and solid lines represent 
forward and bisector fits to the data. The thick 
dashed line shows the trend predicted by the final 
star-formation calibrations of Bell 2003 



on poorly understood physics of cosmic-ray con- 
finement and diffusion operating over 10 7 — 10 8 yr, 
radio continuum surveys have been used to con- 
strain co-moving SFR locally, and at redshifts of 
up to z~l-2 |25I64I84I38| . 

Nevertheless, standard conversion of the radio 
continuum may not properly estimate the SFR for 
all galaxies. Low-luminosity dwarf galaxies, par- 
ticularly BCDs, have higher SFRs than would be 
estimated by applying canonical relations |56| and 
also have globally flatter radio spectral indices, 
clearly reflecting a dominant thermal component 
|5fil2fi| . In these galaxies, particularly those of 
low metal abundance, the SFR derived assuming 
the standard mix of thermal/non-thermal con- 
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tributions is underestimated by a factor of five 
or more |59I5I48| . While rare locally (see [BTJ]'). 
low-metallicity low-mass galaxies might be much 
more frequent at early times and high redshifts, 
given the predictions of the hierarchical merger 
models (e.g. Baugh et al. (this volume), |19|. In- 
deed, such objects may represent the primordial 
"building blocks" - or "sub-galaxies" |2] - in hi- 
erarchical scenarios of galaxy formation. Conse- 
quently, at redshifts z > 3, it is possible that the 
canonical scaling of SFR with radio continuum 
does not apply. Such a systematic mis-estimate 
of SFR could strongly skew our interpretation of 
deep radio continuum surveys. 

3.1. The radio/far-infrared correlation re- 
visited 

While the ubiquity of the radio/far- infrared 
(RFIR) correlation is well documented (e.g. [421 
I24I27I96I70I36I22I74| . it is also clear that the 
slope of the relation depends on the character- 
istics of the sample. Samples dominated by 
low-luminosity galaxies have steeper than unit 
slopes [27170197173] but unit or flatter slopes 
are obtained for high-luminosity samples (e.g. 
|24I22I1()()| ). This is equivalent to saying that the 
ratio of FIR-to-radio luminosity q ( |42] ) : 



Stor Formation Rate 



q = log 



FIR/(3.75 x 10 12 Hz) 



Si. 



4GHz 



(with FIR = 1.26 x 10- 14 [2.58^ 60Aim + ^i o Mm ]) 
increases with decreasing luminosity 3 > as shown 
in Figure [21 

The slope inflection of the FIR-R correlation 
at low luminosities probably stems from a num- 
ber of factors involving both the radio continuum 
and the infrared SED. Moreover, there is some 
evidence that both the radio continuum and FIR 
underestimate the SFR at low luminosities [Hj- 

3.2. Luminosity, metallicity, and age: vari- 
ations in the radio continuum 

As noted above, thermal emission increasingly 
dominates the radio spectrum of galaxies with low 
luminosity and low metallicity in general, and this 
is reflected in their flatter radio spectral indices 
56 55 26 . The same trend is seen in galaxies 
with very young star-forming regions, where non- 




Figure 3. The integrated star-formation history 
of the Universe, from The SFR calculated 

from the "fossil record" in the Sloan Digital Sky 
Survey (SDSS) is shown by the eight large filled 
circles. The other symbols correspond to inde- 
pendent determinations using instantaneous mea- 
surements of the star formation rate at various 
wavelengths. 



thermal radio emission is virtually absent |72|. 
The SFR estimated by assuming a purely ther- 
mal spectrum is roughly ten times higher than 
inferred from a "standard" mainly non-thermal 
spectrum |2*2"] . 

Moreover, optically thick HII regions are very 
likely to be prevalent in young starbursts |59I5I6| . 
This means that instead of a flat thermal spec- 
trum, we would see a rising or highly-absorbed 
spectrum at ^ 1 GHz, which would again alter the 
standard assumptions. Such an effect would also 
be expected to vary with redshift. 

If non-thermal radio emission is suppressed in 
low-luminosity galaxies, the disk magnetic field 
must either be absent or substantially weaker 
than in galaxies with higher luminosities (see 
[56155173) 1. If the starbursts are young, it could 
also mean that there has not been enough time 
to set up the usual 10 s yr diffusion of cosmic ray 
electrons |41|. If the radio emission is also com- 
pact (<200pc), such as that in many BCDs (e.g., 
SBS 0335—052), there can be no ~kpc large-scale 
diffusion. This means that any non-thermal emis- 
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Figure 4. Local (z < 0.1) star-formation density 
for galaxies with star-formation rates between 
0.01 and 1000 M Q yr" 1 (73] . The solid line is de- 
rived from a fit to the local radio luminosity func- 
tion for star-forming galaxies, while filled squares 
show values derived from optical measurements of 
the Ha emission line |34| . In the local universe, 
the most likely place for a star to form appears to 
be in a galaxy with a star-formation rate around 
10 M Q yr -1 . The integral under this plot is the 
zero-point of the 'Madau diagram' which plots 
the integrated star-formation density at different 
redshifts. 



sion must be produced by a different mechanism 
(e.g. evolved supernovae expanding in a dense 
medium |18I2| : young supernova remnants; etc.). 

In any case, the inescapable implication of the 
radio spectral differences as a function of lumi- 
nosity, metallicity, or age is that the conversion 
of radio continuum to SFR may change as a func- 
tion of redshift. While this is probably not an is- 
sue with current facility surveys which probe only 
a limited redshift regime, deep SKA surveys will 
be able to detect not only the most luminous ra- 
dio sources at each redshift, but also adequately 
sample the statistically dominant "normal" star- 
forming population out to redshifts of 10 or so. 
Changes in the conversion factors as a function 
of redshift need to be considered when mapping 
cosmic star formation. 

3.3. The evolving stellar content of galax- 
ies 

The integrated star-formation history of the 
universe (the Madau diagram; see Fig- EJ is likely 
to be accurately mapped out fairly soon. This 
tells us how much star-formation occurred at a 
given cosmic epoch, but to understand galaxy 
evolution we need to know more than this. In 
particular, where and in what range of environ- 
ments, did star formation occur at each epoch? 
Did most stars form in galaxies undergoing "qui- 
escent" star formation, or in galaxies undergoing 
violent starbursts, and how does this vary with 
redshift? How important a role do mergers play 
in triggering star formation? In the local uni- 
verse, star formation occurs in a wide range of en- 
vironments, and in galaxies with star-formation 
rates ranging over at least five orders of magni- 
tude (Fig. . There is no reason to suppose that 
the high-redshift universe is any less diverse. 

Understanding the evolving stellar content of 
galaxies, and the chemical enrichment of the 
interstellar and intergalactic medium which ac- 
companies this, requires a multi- wavelength ap- 
proach. ALMA is sure to have a big impact, and 
will be a powerful tracer of star-forming galaxies 
over a wide redshift range, detecting large num- 
bers of normal star-forming galaxies out to z^3 
and Arp 220-like starbursts to z~10 (e.g. US]). 
But SKA can equally well chart out such lumi- 
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nous star burst galaxies at almost any redshift 
and has the great advantage of a much larger field 
of view and hence superior survey capabilities. 
Moreover, unlike ALMA, the sensitivity of SKA 
to star formation does not depend on the pres- 
ence of metals. Since the first stars must form 
in a primordial metal- free environment, molecu- 
lar clouds would not be expected to contain CO 
or other molecules beyond molecular hydrogen. 
Hence SKA will provide a much " cleaner" less bi- 
ased view of high-redshift star formation. 

4. The assembly of galaxies (0.2 < z < 5). 

SKA has the unique ability to trace the trans- 
formation of HI gas into stars (and into the en- 
riched hot gas which future X-ray missions will 
detect) in a wide range of environments and over 
a large fraction of cosmic time. The fact that 
HI observations trace the depth of a galaxy's po- 
tential, and so can relate star-formation to halo 
mass as well as to luminosity, is also important 
because it allows direct tests of galaxy evolution 
models which deal with mass as well as light. 

4.1. basic SKA capabilities for Hi emission 
observations 

Before discussing how SKA will contribute to 
probing the gas content of the universe over cos- 
mic time we will briefly describe its capabilities 
in terms of sensitivity and resolution to provide a 
baseline for Hi imaging projects. Hi line observa- 
tions typically require a brightness sensitivity of 
a few Kelvin, which in practice limits the angular 
resolution. For SKA this implies the following: if 
one takes as a guideline a Tb sensitivity of ~ 1 
Kelvin in a 0.5" beam at 1.4 GHz (so baselines 
out to at most a few hundred km), then one can 
detect L* galaxies out to a redshift of z — 1.5 
in a 12 hour integration. Although the current 
design goals require a field of view of one square 
degree at 1.4 GHz, a much wider field of view ap- 
pears technically feasible, and would provide an 
impressive gain in survey speed (see also Blake et 
al. this volume). 

Table ^ gives an indication of the 5a Hi de- 
tection limits for a 12 hour integration assuming 
a velocity width of 100 km sec -1 . The calcula- 
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Figure 5. Comoving Hi mass density in neutral 
gas contributed by damped Lya absorbers plot- 
ted as a function of redshift for a cosmology with 
Q M = 0.27, and fl A = 0.73 The region 

fistar is the ±lcr range for the mass density in 
stars in nearby galaxies |33| . The point at z = 
is the value inferred from 21 cm emission from 
local galaxies |l()3j . 



tions assume Ho = 71 km sec -1 Mpc -1 , fiju = 
0.27, and Ha — 0.73. The calculations followed 
the SKA strawman design figures: A rme /T sys = 
20000, 70% of the collecting area is within 100 
km and assumed two polarizations. 

It is clear from Tablc^that it is possible to de- 
tect large galaxies in Hi emission beyond a red- 
shift of z = 1 fairly easily, though with limited 
resolution. For example a galaxy like the Milky 
Way can be detected out to a redshift of about 
1 in a 12 hour integration. Its companions, the 
Magellanic Clouds, can still be seen out to red- 
shifts of a few tenth. Large spiral galaxies such 
as M 101 can be probed much farther, to z « 2.5. 

4.2. The evolving gas content of galaxies 

Various lines of evidence suggest that there is 
strong evolution between z = 1 and z = 3. This 
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Figure 8. Atomic Hydrogen detected in Abell 2192 at z=0.1887. Upper panels: individual channel 
maps from the VLA datacube. The rms noise and Heliocentric velocity are noted above each panel. 
Lower left: Total Hi map. The contours coincide with the position of an inclined barred late- type 
galaxy. Lower middle: contours from the three channel maps plotted on top of each other; solid: 
cz=56483 km/s, dashed: cz=56542 kms, dotted: cz=56600 km/s. Lower right: Ha velocity field of 
the optical counterpart obtained with the PMAS IFU spectrograph at the 3.5m telescope on Calar Alto. 
There is excellent correspondence between the Hi and Ha recession velocities, confirming the Hi detection. 
From Verheijen and Dwarakanath (in preparation). 



is shown not only by the latest estimates of the 
change in comoving star formation density with 
redshift 40 but also by the most recent estimates 
of the evolution in gas density with redshift from 
damped Lya studies |83I37| . 

Direct measurements of the evolution of the 
Hi in galaxies out to redshifts of 1 and beyond 
are not yet available, because present day instru- 
ments lack the sensitivity and resolution to di- 
rectly measure the Hi in galaxies at these red- 
shifts. The largest redshift Hi detection to date is 
a deep VLA observation of the cluster Abell 2192 
at z = 0.1887 |2J which detected Hi in a late 
type galaxy at this redshift (see Fig. 0). Scal- 
ing this observation to the capabilities of SKA 
one expects to get similar results for galaxies at 



distances beyond 2 = 1. 

Damped Lya studies jHS] indicate that the co- 
moving H i mass density is a few times the present 
value beyond z = 1 and out to z = 3 as illustrated 
in Fig. [5] The evolution in comoving Hi mass 
density as shown in Fig. [5] is still dependent on 
assumptions. A recent study of the column den- 
sity distribution function of Hi from observations 
of nearby galaxies by Zwaan et al. [106j suggests 
that there is not as strong an evolution in comov- 
ing Hi mass density as suggested in Fig- El SKA, 
however, will be able to measure the Hi emission 
in galaxies back to redshifts of z ~ 3 directly and 
will revolutionize this area of research. 

There arc various ways in which the SKA will 
contribute significantly to determining the way in 
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Table 1 

Detection limits for Hi emission with SKA W 



z 


Freq. 


Tsys 


Angular ( 3 ) 


Linear 


SB 


Luminosity 


Lookback 


Hi mass ( 4 ) 








Resol. 


Resol. 


Dimming 


Distance 


Time 


limit 




(MHz) 


(K) 


(arcsec) 


(kpc) 


(mag) 


(Gpc) 


(Gyr) 


(M ) 


0.2 


1183.67 


50.4 


0.52 


1.7 


0.796 


0.972 


2.41 


6.1 x 10 s 


0.5 


946.94 


51.4 


0.65 


4.0 


1.486 


2.825 


5.02 


8.7 x 10 8 


1.0 


710.20 


53.8 


0.87 


7.0 


3.026 


6.640 


7.73 


2.7 x 10 9 


1.5 


568.16 


57.5 


1.09 


9.3 


4.000 


11.02 


9.32 


7.2 x 10 9 


2.0 


473.47 


62.7 


1.31 


11.1 


4.796 


15.75 


10.32 


1.5 x 10 10 


2.5 


405.83 


69.6 


1.52 


12.5 


5.469 


20.72 


11.00 


2.6 x 10 10 


3.0 


355.10 


78.3 


1.74 


13.6 


6.052 


25.87 


11.48 


4.3 x 10 10 


3.5 


315.64 


89.3 


1.96 


14.6 


6.566 


31.15 


11.83 


6.7 x 10 10 



1 assuming t = 12 hours, A c /T sys =20000, 2 polarizations and 70% of A c within 100 km) 

2 including a contribution from Galactic foreground emission assuming Tcai(/MHz) = 20( ^° 8 ) 2,7 K 

3 fixed array geometry assumed so that resolution scales with wavelength 

4 assuming 5 rms and 100 km sec -1 profile width. At z=0.2 and z=0.5 the galaxies are assumed resolved 
so here the flux has been added over 8.5 and 1.5 beams respectively. 



which galaxies assemble over cosmic time. Since 
the evolution of comoving Hi gas density as de- 
scribed above has not been directly measured, but 
has been derived from Hi column density mea- 
surements of Damped Lya absorbers it is crucial 
to be able to image the entire Hi disks associated 
with these Damped Lya absorbers. If the associ- 
ated Hi masses are indeed of the order of MJjj = 
6.2 x 10 9 M |104I105I106| then SKA can image 
these disks out to z = 1.4 in integration times 
of 12 hours. Such observations will be crucial for 
linking the existing column density measurements 
to the local Hi density from observations of the 
emission. 

Another way to probe the evolution of galaxies 
over cosmic time is to determine variations in the 
Hi Mass Function as a function of environment 
(local galaxy density) and redshift. From recent 
studies (see Fig. [7] from Zwaan et al. |106| ) it is 
already apparent that the Hi mass function de- 
pends on Hubble type. The Hi mass function also 
appears to depend on environment. The slope of 
the overall HI mass funtion is steeper (~ —1.3 
|1U6| ) than in the Local Group |107| and Ursa 
Major cluster [53] where the slope is ~ —1.0. So 



by measuring the Hi mass function over a range 
of redshifts and in different environments one can 
probe the effects of galaxy evolution. The mod- 
els of Baugh et al. (this volume) also indicate 
changes in the Hi mass function (see Fig.|SJ) which 
SKA will be able to determine. 

In order to do this reliably, sufficient numbers 
of galaxies over a range of Hi masses need to 
be detectable. A deep Hi survey using of order 
30 days of integration time will be able to de- 
tect L* galaxies (which typically have Hi masses 
of 6.2 x 10 9 M assuming H = 71 km/s/Mpc) 
out to redshifts of z w 2.7, i.e. beyond the red- 
shift where the universe shows considerable evo- 
lution. Using the Hi mass function of Zwaan et 
al.(2003) |105| and assuming no evolution of the 
Hi properties of galaxies with redshift one can 
calculate how many galaxies one would expect to 
detect per redshift interval in a one square de- 
gree field of view. Table [2] gives a brief sum- 
mary. These calculations were performed for a 
Qm = 0.27 and S1a = 0.73 cosmology. In such 
long integrations the emission associated with 
Damped Lya systems in the field of view can be 
detected out to z ss 2.7 (assuming they are « Mj^ 
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Table 2 

Detectable Hi Masses as a function of redshift and number of detected galaxies per 
square degree in 360 hour and 1000 hour integrations 



Redshift 


Look Back Time 


Hi Mass Limit W 


Number 


Hi mass limit ^ 


Number 




(Gyr) 




of detections 


(M ) 


of detections 






t = 360 h 




t = 1000 h 




0.5 - 1.0 


5.0 - 7.7 


5.0 x 10 8 


1.6 x 10 B 


3.0 x 10 8 


2.4 x 10 B 


1.0 - 1.5 


7.7-9.3 


1.3 x 10 9 


1.8 x 10 6 


7.9 x 10 8 


2.7 x 10 6 


1.5 - 2.0 


9.3 - 10.3 


2.7 x 10 9 


1.9 x 10 6 


1.6 x 10 9 


2.7 x 10 6 


2.0 - 2.5 


10.3 - 11.0 


4.8 x 10 9 


1.7 x 10 6 


2.9 x 10 9 


2.5 x 10 6 


2.5 - 3.0 


11.0 - 11.5 


7.8 x 10 9 


1.6 x 10 6 


4.7 x 10 9 


2.3 x 10 6 


3.0 - 3.5 


11.5 - 11.8 


1.2 x 10 10 


1.4 x 10 6 


7.3 x 10 9 


2.0 x 10 6 


3.5 - 4.0 


11.8 - 12.1 


1.8 x 10 10 


1.2 x 10 6 


1.1 x 10 10 


1.8 x 10 6 


4.0 - 4.5 


12.1 - 12.3 


2.7 x 10 10 


1.0 x 10 6 


1.6 x 10 10 


1.5 x 10 6 



calculated for the larger z in each interval. Other assumtions as in table 1. 



systems) . 

Out to z ~ 0.3 one will be able to resolve a fair 
fraction of the galaxies so one can obtain rotation 
curves, mass distributions, and gas fractions for 
more than w 10 4 galaxies between now and 3 Gyr 
ago. One then is in a unique position to trace the 
evolution of the ISM in galaxies over a substantial 
fraction of the age of the universe, from the era 
of strongest evolution and star formation activity 
until the present. In addition one would learn 
whether and how the evolution depends on the 
dark matter content and environment. 

The hierarchical galaxy formation simulations 
suggest that most of the evolution in the HI mass 
function is in the mass range Mjji =4 x 10 8 
- 4 x 10 9 M and within the last seven billion 
years. Fig. shows that at z = 1 one expects 
~ 3 times more galaxies in this mass range than 
at z = 0. This figure also illustrates that there 
is very little evolution in the high mass ( Mhi 
> 4x 10 9 M Q ) end of the Hi mass function. This 
is the range where the SKA is especially sensitive 
and can probe easily to at least z = 1, so that 
this prediction can be verified. 

If, on the other hand, the suggestion that the 
gas accretes along filaments and is not heated to 
the virial temperature |10l51l66j is correct, then 
a deep SKA survey should be able to see the fil- 
amentary cool gas reservoirs out to z = 3 since 
until z « 1 there is as much gas in filaments as 



there is in galaxies. Since in this picture the size 
of the filaments that feed the galaxies grow from 
galaxy size at high redshift, to hundreds of kpc 
size structures feeding entire groups at low red- 
shifts, deep Hi surveys should be able to verify 
this idea. If on the other hand mergers of pre- 
dominantly small systems are the main process 
one expects to find a very different distribution 
of Hi structures with redshift. 

In addition, as detailed in the previous section, 
one can use the continuum emission to probe the 
star formation rates of the detected galaxies, in- 
dependent of the effects of extinction using the 
radio continuum-FIR correlation |42I41I21I22I75| 
and test the model predictions. In addition, this 
information can be used to link the star forma- 
tion rates to the Hi contents of galaxies as a func- 
tion of redshift and environment. The global star 
formation rates can be compared with the semi- 
analytic models and the cool gas accretion models 
in order to test the predictions. That environ- 
ment plays an important role is suggested by the 
analysis of Sloan Digital Sky Survey (SDSS) data 
|53I52I45I44| which shows a dependence of galaxy 
mass and star formation rate on environment. 

5. Evolution now (z < 0.2). 

The process of merging and accretion continues 
to the present time. Ample evidence is already 
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Figure 6. The evolution of the Hi mass function 
with redshift in the Baugh et al. (this volume) 
model. The points show the recent determina- 
tion using HIPASS galaxies by Zwaan et al. 106 . 
The curves show the model predictions at differ- 
ent redshifts (solid: z=0 ; dotted: z=l ; dashed: 
z=3, dot-dashed: z=4). 



provided from existing observations, both optical 
and Hi imaging. Zaritsky |l(Jlj and Zaritsky & 
Rix [1(J2| studied several tens of galaxies and as- 
cribed the coincidence of evidence for recent star 
formation from stellar population studies with 
measurable asymmetries in the galaxy shapes, to 
recent accretion events. Kinematic lopsidedness, 
observed in the 2 dimensional Hi velocity fields of 
galaxies |90I85| has also been considered as a re- 
sult of recent minor mergers. Furthermore there 
are at least twenty examples of galaxies which in 
Hi show either signs of interactions and/or have 
small companions |75l76l87j . This suggests that 
galaxies often reside in an environment where ma- 
terial for accretion is available. An example is 
given in Fig. [5] which shows a new Hi image of 
the galaxy NGC 3359 which is accreting gas from 
a faint Hi companion. 

Other clear evidence for ongoing minor merger 



Figure 7. The Hi mass function for different mor- 
phological types. Shown are Schechter fits (solid 
lines) and the overall Hi mass function for the to- 
tal sample (dashed line). From Zwaan et al. |106| . 



activity comes from the discovery of the Sagittar- 
ius dwarf galaxy [4T)| and of substructure in the 
halo of the Local Group galaxy M 31 [50123163] . 

The Hi Rogues Gallery of Galaxies |43| . 
www.nrao.edu/astrores/HIrogues/) gives a 
glimpse of the role of gas in galaxy evolution 
at z = 0. Ideally one would want to make com- 
plementary Rogues Galleries at z = 0.5, 2 = 1, 
z = 1.5, and z = 2 to study the process of accre- 
tion in detail. Inspection of the Rogues Gallery 
clearly demonstrates the importance of using the 
Hi for probing galaxy evolution. The prospect 
of obtaining similar information to either higher 
depth in sensitivity or to much larger redshifts is 
truly exciting as radio observations are the only 
means to probe these aspects of galaxy evolution. 
There simply is no other way to trace the neutral 
atomic gas around galaxies. A very good illus- 
tration is the extended Hi around SO galaxies 
in the Ursa Major cluster shown in Fig. ^| SO 
galaxies are thought to be the result of galaxy 
transformation by gas removal (stripping and/or 
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^GC -till 




Figure 10. Extended Hi filaments observed near the brightest lenticular galaxies in the Ursa Major 
cluster. These lenticulars are located in small sub-groups within the Ursa Major volume, and the 
Hi morphology indicates that tidal interaction between the galaxies is taking place. From Verheijen and 
Zwaan [§T] 



tidal effects) but in these objects this is not so ob- 
vious. A second example of the diagnostic power 
of Hi imaging is NGC 4522 in the Virgo cluster, 
where the Hi morphology presents overwhelming 
evidence that stripping by ram pressure is taking 
place, (see Fig. ITO 

That indeed the environment plays an impor- 
tant role in determining the Hi properties (and 
hence probably the evolution) of galaxies has be- 
come clear from different lines of evidence (for a 
recent review see van Gorkom |89j): the finding 
that galaxies in clusters out to 2 Abell radii are 
deficient in Hi as compared to galaxies of simi- 
lar type in the field (SHI) see also Fig. ^| ) and 
the truncated and sometimes displaced Hi distri- 
butions of galaxies in rich environments (result- 
ing from ram pressure stripping and tidal effects). 
Additional evidence for strong environmental ef- 
fects is shown in Fig. 1131 from Verheijen et al. 
|91| . In the denser environments of the core of the 



Virgo cluster the Hi deficiencies and Hi to opti- 
cal sizes are notably different (respectively larger 
and smaller) than in the less dense environment 
of the Ursa Major cluster and the field. 

The velocity structure and spatial distribution 
of Hi detected galaxies near clusters suggests 
that infall is still going on and that clusters are 
still capturing galaxies which fall in along large 
scale structure filaments. Evidence such as the 
Hi structures around SO galaxies and the stripped 
gas in NGC 4522 (Fig. Q33 and Fig. Q3J are cru- 
cial in pinning down these processes and SKA will 
be able to see such structure over a considerable 
range of redshift (i.e. evolution). 

A shallow survey covering 100 square degrees 
in two months (assuming 100 km sec -1 profile 
width, a 5<T detection limit, a field of view of 
one square degree at 1.4 GHz , and an instan- 



1 if the SKA has a field of view of 50 square degrees rather 
than the one square degree of the SKA strawman design, 
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NGC 3368 HI overlay on DSS 
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Figure 9. Hi distribution of NGC 3359 at a reso- 
lution of 30" superposed on the digital sky survey 
image (left panel). Contours are 0.1, 0.2, 0.4, 0.8, 
1.6, 3.0, and 5.0 xlO 21 cm- 2 [H?| 
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Figure 11. NGC 4522 in the Virgo cluster from 
Kenney et al. 5 1 . The Hi contours are displaced 
from the optical disk and clearly show the sweep- 
ing effect of ram-pressure stripping. 



taneous bandwidth of 500 MHz) will be able to 
see NGC 3359-like structures to z = 0.15 opening 
the exciting prospect of being able to probe the 
environment of several million galaxies and pin 
down the merger/accretion characteristics at the 
present time over a range of galaxy densities. The 
resolution and sensitivity of such a survey will be 
comparable to that of the images shown in Fig. 
and Fig. 1101 The volume covered by such a survey 
is ~ 15,000 Mpc 3 and will detect ~ half a million 
galaxies, providing a detailed Hi image, not only 
of the galaxies but also of their Hi environments. 
Only with such data can the evolution of galax- 
ies, in particular the mechanisms for converting 
gas into stars in and around galaxies be studied. 
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